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Abstract 

We show that any scalar differential operator with a family of polyno- 
mials as its common eigenfunctions leads canonically to a matrix differen- 
tial operator with the same property. The construction of the correspond- 
ing family of matrix valued polynomials has been studied in |D1UD2IIDV] 
but the existence of a differential operator having them as common eigen- 
functions had not been considered This correspondence goes only one way 
and most matrix valued situations do not arise in this feishion. 

We illustrate this general construction with a few examples. In the 
case of some families of scalar valued polynomials introduced in |GH) we 
take a first look at the algebra of all matrix differential operators that 
share these common eigenfunctions and uncover a number of phenomena 
that are new to the matrix valued case. 



1 Introduction 



The last few years have witnessed some progress in the problem of finding ex- 
plicit families of matrix valued orthogonal polynomials that are joint eigenfunc- 
tions of some fixed differential operator with matrix coefficients, and a (possibly 
matrix valued) eigenvalue that depends on the degree of the polynomial. Most 
of these examples involve differential operators of order two, an issue raised in 
the case of symmetric operators in [D3 . The subject of matrix valued orthog- 
onal polynomials, without any reference to differential operators, was initiated 
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in two papers by M.G. Krein, |K1[ IK2| . The search of situations where these 
polynomials satisfy extra properties, such as the one above, makes it more likely 
that they will be useful in applications. 

There are by now four main ways to search for cases with this particular ex- 
tra property: appealing to some group representation structure, as in [GPTli 
IGPT2[ [GPT31 IGPT4[ [FT) . solving (in the case of sy mmetric diff erentia l opera- 
tors) an appropriate set of differential equations as in [DGl[[DG2llD4[[GllGPT3j 
or a set of moment equations as in [D3j . and finally and so far much less suc- 
cesfully, solving the so called ad- conditions as in |CG3| . One can see that these 
ad-conditions are necessary and sufficient for the existence of a differential op- 
erator, but finding all solutions remains a huge task. Most of the examples 
obtained by this route correspond most likely to an "orthogonality functional" 
that is not given by a positive definite matrix as envisaged by M. G. Krein. It 
appears reasonable at this point to consider in some detail any method that 
gives explicit examples. 

One such path is taken in [CGlj where one finds instances where there is a 
differential operator of order one, and in |CG2j where one explores the algebra of 
all differential operators going along with a fixed family of differential operators. 

In this paper we take the initial steps in yet another path. This comes about 
by putting together two previously unrelated lines of work. 
The first one comes from studying the relation between scalar valued polyno- 
mials satisfying higher order recursion relations and matrix valued polynomials. 
This has been considered in [Dl|, ID2( IDVj . and is reviewed here in section 2. 
These papers do not look into the issue of possible differential equations satisfied 
by these polynomials. 

The second one comes from an extension of the well known Krall-Laguerre 
and Krall-Jacobi orthogonal polynomials. This has been considered in [GHj 
to produce, via Darboux factorization, families of (scalar) polynomials that 
(in general) satisfy five term recursion relations and fourth order differential 
operators with polynomial coefficients. The relevant parts of this construction 
are recalled in section 4. 

These two tools are combined in this paper as follows: we first establish a very 
general result that produces one matrix valued differential operator every time 
that we start with a scalar one (but not necessarily the other way around, 
as spelled out later). This general result is first illustrated in the case of the 
Laguerre polynomials where the orders of all operators is equal to two. This 
completes the material in section 3. We then illustrate, in section 4, the case of 
higher order operators by looking at the extended scalar valued Krall-Laguerre 
and Krall-Jacobi polynomials. 

A sligthly more detailed description of the contents of sections 3 and 4 follows. 
Given specific families of scalar valued polynomials we consider the correspond- 
ing family of matrix valued polynomials as in [DH ID2| IDVj . We find that in the 
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Krall-Laguerre case there is (up to scalars) only one differential operator with 
matrix valued coefficients (and of order four) going along with each one of these 
examples. In particular, there is no lower order operator with this property. By 
going up to order six we see that the corresponding algebra cannot be generated 
by the fourth order operator constructed in the canonical fashion considered in 
the first part of the paper. The case of the extended Krall-Jacobi polynomials 
leads to a richer variety of examples. In particular we get a third order dif- 
ferential operator in the matrix case which cannot yield a third order scalar 
differential operator. All of this is discussed in several subsections of section 4. 

To better compare the results of this paper with previously known material 
notice that the situations uncovered so far feature an algebra that includes dif- 
ferential operators of low order, typically one or two. There are even situations 
where one finds a nontrivial (i.e. not a scalar multiple of the identity) differential 
operator of order zero. One could wonder as to the existence of situations where 
the lowest order operators in the algebra needs to be of order higher than two. 
The group theoretical line of attack mentioned above, is not likely to produce 
examples of this kind, since in this case there is a Laplace-Beltrami operator 
around. Solving the appropriate differential equations for the higher order case 
appears to be a difficult task, while solving the corresponding ad-conditions 
appears entirely hopeless at this point. 

As a by-product of the results in this paper we produce examples of the type 
envisaged above, which as far as we know have not been seen before. 

Another new result is the observation that we have here a one way street; while 
the results in [Dl|, ID2| IDV| establish a reversible link between scalar valued and 
matrix valued polynomials we see here that if the scalar family is made up of 
common eigenfunctions of some differential operators this property extends to 
the matrix valued polynomials but not (necessarily) the other way round. 

We close this introduction with some remarks, aimed at a broad audience, about 
some of the possible uses of matrix valued orthogonal polynomials. 
One of the many applications of the theory of scalar orthogonal polynomials is 
to a detailed study of a special kind of Markov chains with state space given 
by the non-negative integers and with a one step transition probability matrix 
that allows only moves to the nearest neighbours. The basic paper here is 
|KMcG| . These authors derive an expression (based on the spectral analysis 
of the corresponding difference operator) for the n-step transition probability 
between states i and j in terms of the orthogonal polynomials going along with 
the one step transition probability matrix. There are several advantages to such 
a formula. For instance if one knows the orthogonality measure one can compute 
the probability of going from states i,j with i < j in any number of steps from 
the knowledge of the first j rows of the one step transition matrix. Other uses 
of this connection with orthogonal polynomials lie a bit deeper: by using the 
Stieltjes transform of the orthogonality measure one can see that (in principle) 
the knowledge of the probability of going from state to state in an arbitrary 
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number of steps determines the full n-stcp transition probability matrix for any 
n. 

In the last section of the paper by Karlin and McGregor mentioned above, these 
authors consider the problem of using their methods , based on a three term 
recursion relation that is semi-infinite, to the case of a doubly infinite three 
term recursion. Although they stop short of introducing explicitly the notion of 
matrix valued orthogonal polynomials they compute in detail the orthogonality 
matrix measure that is relevant to this problem. The full blown analysis of the 
problem of a random walk (with nearest neighbour transitions) on the set of all 
integers that was the concern of the last section in [KMcG| has been done in 
detail later in [DRSZj . |G1| . This extension from the study by means of scalar 
valued orthogonal polynomials of birth-and-death processes to the study of so 
called Quasi-birth-and-death processes by means of matrix valued orthogonal 
polynomials is already a concrete example of a good motivation for this study. 
For a situation where the matrix valued orthogonal polynomials happen to sat- 
isfy matrix valued differential equations which along with the orthogonality 
measure are an outgrowth of the work started in |GPT1| , see |GIj . 

The best known families of scalar valued orthogonal polynomials, connected 
with the names of Hermite, Laguerre, and Jacobi owe their importance to the 
differential equations they satisfy and thus they appear in connection with sim- 
ple situations involving the Schrodinger equation. In the scalar case the only 
possible examples of orthogonal polynomials with a differential operator float- 
ing around are the ones mentioned above. On the other hand the situation in 
the matrix valued case is much more complicated and it opens the door to an 
embarrassment of riches in terms of examples. The study carried out here is an 
effort to produce more examples of this novel situation as a step towards other 
applications. 

2 Scalar and matrix valued polynomials 

The relation between scalar polynomials satisfying higher order recursions and 
matrix valued ones satisfying three term recursions was found by one of us (see 
[m] and [D2]) and it is well described in [DV] . 

The procedure shows a way to produce a vector polynomial from an scalar one 
p{x) by considering the decomposition of p{x) congruent with the residues of 
the exponent of x in p{x) modulo N. 

This is done by using the operators RN.m, m = 0, • • • , iV — 1, defined by 



i.e., the operator Rn,™. takes from p those powers congruent with m modulo N , 
then removes the common factor 2;™ and changes x^ to x. 



(2.1) 
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From a sequence of scalar polynomials {pn)n we can then define a sequence of 
matrix valued polynomials Pn{x), n — 0,1, ■ ■ ■ , by the recipe 



(2.2) Pnix) = 



( RN,o{PnN)ix) 
RN,oiPnN+l)ix) 



RN,N-liPnN)ix) ^ 
RN,N-liPnN+l)ix) 



\RN.,o{PnN+N-l){x) ... RN.N-l{PnN+N-l){x) J 

This process can be reversed since for any polynomial p{x) we have 

p{x) = i?w,o(p)(a;^) + xRN,iip){x^) + ... + x^-' RN^N-iip^x"") 

It turns out that if Pn(x), n ~ 0,1,2, ... , satisfy the {2N + l)-term recurrence 
relation 



N 



(2.3) 



X^Pn{x) = Cn,OPn{x) + ^ [c„, _fep„-fc (a;) + Cn,kPn+k{x)], 



k=l 



{pn{x) vanishes if n is negative), then the matrix valued ones P„, n = 0, 1, 2, . . . , 
satisfy the following three term matrix recurrence relation 



(2.4) 



xPn{x) = AnPn+i{x) + BnPn{x) + C„P„_i(a;) 



where the matrix coefficients are nothing but the N x N blocks of the (2iV+ 1) 
banded matrix associated to the scalar family (p„)„ featuring in the higher order 
recursion that defines the Pn{x). This semi- infinite banded matrix has rows that 
look as follows 



first row 



second row 



third 



C0,0, Co,l, • • • , Co,Af, 0, 0, 0, ■ 



Cl,-1, Ci,o, Ci 1 , Ci 2, ■ • • , CiJV, 0, 0, 0, • 



C2,-2, C2,-l, C2^0, C2,l, C2,27 C2,3i ' ' ' j C2,Af 



and for i large enough, the {i + 1)^^ row looks like 

0, 0, . . . ,0, Ci,_JV, Ci,-Af+1, ■ ■ • , Ci^-i,Cifl, Ci,i, . . . , Ci,N,0, 0, . . . 

Given a block-tridiagonal matrix (that defines, by using (|2.4p . a sequence of 
matrix valued polynomials) one can always assume, by using an n dependent 
unitary matrix, that the off-diagonal blocks are triangular, and thus the matrix 
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becomes a banded matrix with scalar entries. This shows that using the opera- 
tors RN,m, TO = 0, • • • , -/V — 1, one can also go from matrix valued polynomials 
satisfying a three term recurrence relation like (|2.4p to scalar valued polynomials 
satisfying a higher recurrence relation as in (|2.3p . 

The polynomial in the formula (|2.3[) does not play any special role. One can 
replace it by (x — a)^ for an appropriate choice of the constant a; in this case, 
one has to use the operators 



instead of the RN,m defined in (|2.ip (see |DV| . pp. 269, for the general case 
when in (|2.3p . the polynomial is substituted by a generic polynomial h{x)). 

We say that a higher order recursion (|2.3p for {pn)n is symmetric if c„._fc = 
Cn-k,k, and, analogously, that the matrix three term recurrence relation (|2.4[) is 
symmetric if = C*_^_i. As it was shown in [DVj . a symmetric higher order 
recursion gives a symmetric matrix three term recurrence relation. 

The symmetric case corresponds to a Hermitian (2iV + l)-banded matrix and 
the associated matrix polynomials defined by (|2.2|) are then orthonormal with 
respect to a weight matrix (i.e. a positive definite matrix of measures having 
finite moments of any order). 

Given a higher order scalar recursion that is not symmetric there are two ways 
to try to symmetrize it. Starting from the polynomials Pn{x) defined by the 
recursion one can define new polynomials r„(a;) — ^nPn{x), and look for a 
convenient choice of the sequence of numbers (7„)ti so that the corresponding 
higher order recursion for (r„)„ is symmetric. 

A different, and more general approach consists of producing from the scalar 
polynomials Pn{x) the matrix polynomials Pn{x) as above. Now one consid- 
ers a new sequence of matrix polynomials Rn{x) — TnPn{x), and looks for a 
convenient choice of the sequence of matrices (r„)„ so that the matrix three 
term recurrence relation for (i?„)„ is symmetric. It is easy to see that the first 
attempt is a special case of the second one, which consists of insisting that the 
matrices r„ be diagonal. 

3 Producing matrix differential equations out of 
scalar ones 

Our starting point is a sequence {pn)n of polynomials, with deg p„ = n, satis- 
fying the differential equation 




m 



(3.1) 
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where ai, I — 1, • • • ,m, are polynomials of degree not bigger than I and 7„ are 
real numbers. 

We split up the polynomials Pn as explained in the previous section to get a 
sequence of matrix polynomials P„ of size N x N. We now show that these 
matrix polynomials inherit a matrix differential equation from the differential 
equation (|3.ip for the scalar polynomials p„. The proof is constructive, so that it 
produces an explicit expression for the differential coefficients Ak, k — 0, - ■ ■ , m, 
in the new differential operator. 

For the sake of simplicity, we show the result in full only for the case N — 2 
(the general case can be proved in the same way). 

Theorem 3.1. Assume that (p„)n is a sequence of scalar polynomials satisfying 
a differential equation of the form i3.1\) where ai, I — 1, ■ ■ ■ ,m, are polynomials 
of degree not bigger than I and 7„ are numbers. Then the matrix polynomials 
Pn defined from pn by satisfy the differential equation 



fc=0 

where the coefficients A^, k — 0, - ■ ■ ,m, are given by 



Ak = (^Cfe 
Ck,iix) 



l=k J 



bk,Loi^) bk,Li{x) \ f ais){x) ai^i{x) 

xbk,i,i{x) +lbk.i-ifi{x) bk,i,a{x) + lbk,i-i.iix) J \xais{x) aifi{x) 



bk,i^ E {-iy(%j{2j-l)---{2j-l + l)\x 



2k-l 



(as usual [x] denotes the integer part of x) and we write po — R2,o{p) o-'^^d 
pi = i?2,i(p), where i?2,o o.'fT'd R2.1 o,i"e the operators defined in h2. 1\] for N — 2. 
The eigenvalues r„, n >0, are given by 



(inN 
7nAf+l 



V 







InN+N-lJ 



Proof. To simplify the expressions, we again write pq — i?2,o(p) and pi 
R2,iip), so that p{x) = poix'^) +xpi{x'^). 
A simple computation shows that 



fc=[(i + l)/2] 
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where 



bkA^)=[^-^ J2 {-iyr.)2j{2j~l)---{2j-l + l)]x''^-\ k>l/2. 

j = [(/+l)/2] ^-^^ 



Applying this formula to Pn{x) — Pn.o{x'^) + a;p„,i(x^), we get after some sim- 
plifications that 

I 

= J2 {bkM^')pi%{x') + {xHk,iA^') + ibkj-M^')) pi%')) 

\fc=0 / 

where we take bk^k' = 0, when k > k' . This can be written in the form: 
(3.2) (p„,o(x2)+a;p„,i(x2))W 



k=0 



ybk,Li{x^) + lbk,i-i,oix^)) 



If we write 



bk.i,i{x) \ f ai^x) ai^x) 

+ lbk,i-i,i{x) J \xaLi{x) ai^x) 



C (x) = ( bk,ifl{x) b 

''•^^ \xbk,i,i{x) + lbkd-i,o{x) bk,ifl{x) 

and take into account that for any polynomials q one has {pq)o = paqa + xpiqi 
and {pq)i — Poqi +Piqo, we get from p.2p that 

(3.3) [[ai{x)pli\x))^ (a,(x)p«(x)) J = ^ (^^(x) p^^lix)) C^,. 

k=0 

From here, it follows that the differential equation (|3.ip gives the following 
matrix equation for the vector (pl^l{x) p^^\{x)j : 



fc=0 \l=k 

It is now easy to finish the proof. 



□ 
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We illustrate the general construction above with an example. 

Take the Laguerre polynomials {L")n, which satisfy the second order differential 

equation 

x{LZ{x))" + ia + l- x)iLZix)y = -nLZix). 
The proof of Theorem 13.11 gives for the matrix polynomials 

2 2y/x 

L2n+liV^ + -^'2n+l(~\/^) -^2ji+l (v^) ~ -^2n+l (^V^) 



\ 



2Vi 



the matrix differential equation 



(-^n,2x2) 



Ax 
4x2 Q 



-2n 




-2x 2a + 4 
(8 + 2a)x -2x 



l: 



n,2x2 





a + 1 -1 





'2n - 1 



TO 



In this case, the sequence (^",2x2)" is orthogonal with respect to the weight 
matrix 

■ 1 ^ 

\/x 

1 v^y 



For any size iV, the matrix polynomials (L" nxn)" defined from the Laguerre 
polynomials by (|2.2[) satisfy the matrix second order differential equation 



(Knxn)" Mx) + {LIn>cn)' M^) + Ll2x2Ao 



(inN 
7riAr+l 

V 



\ 





InN+N-lJ 



L 



n,2x2: 



where 



Aiix) 









( 





• • 

















4x2 





• • 





















4x2 


• • 














A2{X) 






























• • 


















I 





• • 


• 4x 








( 


— X 



















4 + a/2\ 


(iVH 


h 2 + a)x 




— X 






















(iV + 44 


- a)x 


— X 































— X 




V 

















{N + 


2(A^-1)H 


- a)x —X J 
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and 



Aoix) 



(N-l/N'^ 

a + l/N^ N-2/N^ 

2a + N -i/N^ 



V 



{N- 







l/iV2 

l)a+ (A^- 



1)V7V2 








0/ 



4 A look at some examples involving higher or- 
der operators 

In section 2 we have shown that a scalar valued sequence of polynomials that are 
common eigenfunctions of a given differential operator gives rise, in a canonical 
fashion, to a sequence of matrix valued polynomials and to a differential operator 
that has them as common eigenfunctions. 

The purpose of this section is to illustrate a number of issues suggested by this 
general result in the case of some examples. These examples deal with what 
are called extended Krall-Laguerre and extended Krall-Jacobi polynomials in 
[GH]. In both cases the scalar polynomials we consider satisfy a unique fourth 
order differential equation which therefore gives rise to a fourth order differential 
operator with matrix coefficients. The general theme of this section centers 
around the existence of other matrix valued differential operators having these 
same polynomials as common eigenfunctions very much in the spirit of [CG2] . 

4.1 The scalar extended Krall-Laguerre polynomials 

The scalar valued (extended) Krall-Laguerre polynomials first introduced 

in [GHj are given, for n = 0, 1, 2, . . . , by 



where 
(4.1) 
(4.2) 



Pn{x) = -{qn+l(x) +Xn+iqn(x) + J/„g„- 1 (x) ) 
X 



2n^ + 2{R - l)n - R 



Vn = 



{n -1) + R 
n{n + a){n + 1 + i?) 



n + R 

and qn{x) = i"(a;) stand for the familiar Laguerre polynomials given by 
L$J(x) - {x - (2n - 1 + a))Ll_^{x) - {n - l){n - 1 + a)Ll_^{x) 
and L-^ii^x) = 0, L^{x) = 1. 

Here i? is a free parameter, and so is a, which is taken to be larger than —1. 
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The Pn{x) are obtained in [GH] by an application of the Darboux process start- 
ing from the Laguerre polynomials which are orthogonal and satisfy the 
three term recursion relation given above. For generic values of a > —1 the 
Pnix) do not satisfy such a three term recursion relation, but they still satisfy 
the five term recursion relation given by 



(4.3) 
fxi 1 

yi X2 1 

2/2 a;3 



fxi 1 

yi X2 1 

m X3 



\ 



fpo{x)\ 
Pi{x) 

P2{x) 



fpo{x)\ 

Pi{x) 

P2{x) 



V 



7 



V 



7 



V • / 



V • / 



where x„ and t/„ are given above and x„ and are given as follows 

2n2 + 2{R - l)n - R 



n + R 
n(n + a){n — 1 + R) 



(4.4) x„ 

(4.5) yn 

^ ' ^ n+R 

The point of the construction in |GH| is that the new family of polynomials 
Pnix) are common eigenfunctions of a fixed fourth order differential operator in 
the spectral parameter will be seen shortly. 

The method discussed in |GH| shows that a proper use of the Darboux process 
will take us from a trivial bispectral situation, i.e. the Laguerre polynomials 
L^{x) are eigenfunctions of a fixed second order differential operator and they 
are the common eigenvectors of a tridiagonal difference operator, to a much less 
trivial one: the polynomials Pn{x) feature in a higher order bispectral situation 
that is not a trivial consequence of the original one. 
In our case we have 

(4.6) 

[x^D'^ - 2x{x -a- 2)D^ + {x{x - 2{R + a) - 6) + a{a + 3))!?^ 

+ 2{{R + l)x -{a + 1)R - a)D + R{R + l)]pn{x) 
= {R + n){R + n+l)pn{x). 



We now produce out of these Pn{x) a sequence of matrix valued polynomials of 

size 2x2 using the method explained in Section 2. 

For concreteness we display below Pq{x) and Pi{x). We have 



Po{x) 



Pi{x) 



1 

a+{a+l)B. 
R+1 



(a+2)({a+l}R+2a) 



—3x 



((a+3)fl+3a+8 
R+3 



R+2 

(a+2)(a+3)((a+l)J?.+3a) 
R+3 



2({a+2)R+2a+3) 
R+2 

3(a+3)((a+2)Ji;+3a+4) 
R+3 



Here we display a fourth order differential operator B, arising by applying the 
method explained in the previous section, that has the matrix valued polyno- 
mials {Pnix))n as a set of common eigenfunctions 
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B = D^x^I + D 



D 



— X 

[a + 7)a;^ 

(2J?.+3)a _|_ a(ct+3) 
^ fa+4)fl+4(a+3) ^ 



D 



j a;^+3:(a +9a + 15) {R+a+6)x 
2/4 2 

_ CR+£+9)£^ a:^ + (Q^ + 15a+39)3: 
2 4 , 



(a+2)fl+2a+3 

a+4)fl+4(a+3) 2fl+5 , 3(^0 + 1) (a+4) 
4 8 8 , 



__R+1 

(ct+l^fl+g 







The matrix valued eigenvalue is read-ofF from the expression 



Pn{x)B 



R+l 



o\ , 



0; ' V 



2_R.+3 
8 



n+—I 
4 



Pnix) 



Notice that our differential operator acts on the right on its argument. This is 
also the reason why the eigenvalue matrix appears before the eigenfunction P„. 

Observe that the differential operator above depends both on R and a while 
the eigenvalue does not. If one were to consider, for a fixed pair R,a the 
family of matrix valued polynomials, and then in the spirit of |CG2| the full 
algebra we would end up getting isomorphic algebras of differential operators 
with matrix coefhcients for different values of a. This follows from the fact that 
these algebras of differential operators are isomorphic to the algebra of matrix 
valued eigenvalues. 



4.2 Matrix extended Krall-Laguerre and orthogonality 

The sequence of matrix polynomials (P,i)„ is not always orthogonal with respect 
to a weight matrix of measures. This only happens if each polynomial P„ can 
be normalized using a nonsingular matrix: Qn = r„P„ so that {Qn)n is an 
orthonormal sequence with respect to a weight matrix. In that case, the sequence 
{Qn)n satisfies a three term recurrence relation of the form 

xQn — ^n+lQn+l + B„Qn + A^Qii-l, 

with An non singular and _B„ Hermitian. 

If we were dealing with a monic sequence of polynomials P„ satisfying a three 
term recursion, such as 

xPn = Pn+l + BnPn + ^nPi-l, 

this symmetrization can be done only when one can find positive definite ma- 
trices Sn such that -B„S'„ is Hermitian and AnSn+i = Sn- 

As we noticed at the end of section 2, given a scalar family of polynomials 
Pn (x) we could attempt a scalar symmetrization of the (higher order) recursion 
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relation. Given the Krall-Laguerre polynomials (p„)„ we can look for a normal- 
ization r„(x) = TnPn{x) satisfying a symmetric five term recurrence relation of 
the form 

where a„ 0, n > 0. 

If we use the parameters (a;„)„, (y„)„, (x„)„ and (?;„)„ (see (|4?T|) . (|4?2]) . 
and (|4.5p ) it turns out that we would need 

_ _ {VnXn + Xn+iyn){yn+lXn+l + Xn+2yn+l) 

This only happens when a = and then: 



a, 



.+2,i?,o = (n+ l)(n + 2) 



(n + i?)(n + i? + 3) 



6n+i,-R,o = 2(n + 1) 



^ (n + i?+l)(n + i? + 2)' 
2(n + + (2n + l)(2n + 3)i? + 2n(n + l)(n + 2) 



(n + i? + l)^(n + R){n + R + 2) 

_ n'^{n + R+lf + {n + lf{n + Rf + ((2n(n + 1) + (2n + 

~ (n + i?)(n + i?+l) ■ 

For a = 0, the five terms recurrence relation (14. 7|) is just an iteration of the three 
term recurrence relation for the Laguerre type orthogonal polynomials {L^)n 
introduced by Koornwinder (see [Kooj and also |KoKo| ) . They are orthogonal 
in [0, +00) with respect to the measure + ^Sq. The four order differential 
equation (|4.6p reduces to the well-known four order differential equation for this 
family (see, for instance, |LKj . p. 112). 

The corresponding matrix polynomials {Pn,B.,o)n are then orthogonal in [0, +00) 
with respect to the weight matrix 




In this case one can see that the more general attempt to symmetrize the matrix 
valued recursion relation fails too, except in the case a = 0. 



4.3 The scalar extended Krall-Jacobi case 

The scalar valued (extended) Krall-Jacobi polynomials Pn{x), first introduced 
in [GH] are given by 

1n{x) = ]—-{ynPn-l{x) + Xn+lPn{x) + Pn+i{x)) 

X + 1 
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where, as in (4.11) and (4.12) in [GH], wc have 

_2{a + (3 + n){aP + (3^ + 2na + 2{n - l)(3 + 2n{n - 1)) 
^" ~ (a + /? + n - l){a + f3 + 2n - 2){a + f3 + 2n) 

2{a + n- 1)R 
^ {a + P + n-l){a + P + 2n- 2)6'„_i 

On+l 



and 

9n = + (a + f3)n + R. 
Here Pnix) stands for the usual (monic) Jacobi polynomials satisfying 

XPnix) = a„Pn-lix) + bn+lPn{x) + p„+i{x) 



with 



4n(n + a + l3){n + a){n + (3) 

(2n + a 

&n = 



(2n + a + /3 - l)(2n + a + /5)2(2n + a + /? + 1) 



(2n + a + /3 - 2)(2n + a + /5) ' 

As can be seen in [GH^ , the crux of the construction of these polynomials is that 
the polynomials p„ (x) and ix) are also related by means of 

V-a{x) = — \—{ynqn-\(x) + S„+i(7„(a;) + g„+i(a;)) 
a; + 1 

As explained in [GH] the expressions for a;„ and j/„ can be obtained by replacing 
hn by hn + 1 in (3.17) of that paper. The resulting expression for Xn is correctly 
given by the first part of (4.13), but there is an error in the expression for y„ 
reported in [GHj . and we give the correct expression below, namely 

Vn = a;„+i(x„+2 - hn+\ - 6n+2 " 2) - ?/„+i + a„ + a„+i + (6„+i + Xf 
For completness we reproduce the expression for a;„ from [GH] , namely 

_ 2(a + /3 + n - l)(a/3 + /^^ + 2(n - \)a + 2n/3 + 2n(n - 1)) 



+ 



(a + /3 + n)(Q; + /3 + 2n - 2)(q; + /3 + 2n) 
2(a + n)i? 



(Q; + /? + ?i)(a + /3 + 2n)6'„ 
Finally i? is a free parameter. 

These polynomials qn{x) satisfy, for generic values of a, (3, a five term recursion 
relation and they are the common eigenfunctions of a fourth order differential 
operator PQ with 

P = (a;2 - 1)1)2 ^{^{^a + P+Vjx- l3 + a + l)D + R 
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and 



Q = (a;2 - 1)D^ + {{a + p + 3)x - p + a - 1)D + {R + a + p + 1). 



More explicitly we have 



with 



Xn = {R + n{a + P + n)){R + (n + l)(a + /3 + n + 1)). 



The five term recursion relation expresses the product 

as a linear combination of the polynomials qn+2, Qn+i,Qn, Qn-i,Qn-2- 

One could raise again the possibility that this five recursion could be made 
symmetric by premultiplying the scalar polynomials by appropriate n depen- 
dent constants. This would require, as in the case of the modified Krall- 
Laguerre polynomials that a certain identity should hold among the quantities 
Xn,yn,Xn,yn- An explicit computation using the expressions above shows that 
this happens only for /? = 0. 

If one obtains the matrix valued polynomials out of these scalar ones and tries 
to symmetrize the resulting three term recursion relation one finds, once again, 
that in this case nothing is gained by allowing this more general cure to our 
problem: there is a positive definite orthogonality matrix only in the case of 



We consider now (as an example) the special case given by a = 3/2, /3 = 7/8, 
R = 7. In this case the scalar fourth order operator mentioned earlier becomes 



The general method indicated earlier produces a fourth order differential opera- 
tor with matrix coefficients satisfied by the matrix valued polynomials obtained 
from the sequence Pnix). In the spirit of [CG2] one can study the algebra of all 
differential operators associated with this family of matrix valued polynomials. 
In this example we have checked that the algebra contains only one (up to 
scalars) differential operator of order four. If we look for operators of order less 
than or equal to six, we find only one such operator (modulo operators of order 
4), whose highest order terms are 



/3 = 0. 



(a;2 - l)(51a;-|-5) 



7 35 
D^ + — {503x^ + 98a; - 281)£)^ + — (83a; + 21)D. 



4 




3 5 /(a; - l)2.T2(107a;-40) 
W V 5{x-lfx^ 



5(a;-l)2a;2 
(a; - 1)2x2 (123a;- 56) 



+ ... 
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In particular the algebra cannot be generated by the fourth order operator 
constructed out of the scalar one. As mentioned in the introduction, this also 
gives examples of situations where the lowest order of a (nontrivial) differential 

operator in the algebra is higher than two. 

Extensive computations suggest that these two observations hold for generic 
values of a, 0, R. Just as in a previous example, and all the examples below, 
all differential operators with matrix valued coefficients act on the right of their 
matrix valued arguments. 



4.4 Lower order operators connected to the extended Krall- 
Jacobi case 

A remarkable phenomenon takes place when a and /3 differ by one. Consider the 
case /3 = a + 1. Here, very much as in the generic case, we find that the scalar 
extended Krall-Jacobi polynomials are the common eigenfunctions of only one 
fourth order differential operator, namely 

(a;2 - ifD^ + 2(1 - x'^){{2a + 5)x - 1)D^ 

- 2(a;((2a^ + 11a + 19)a; - 2a - 2) - 3a - 13)Z)^ 
- 2((2a + 3)(2a + 9)a;- 7)£). 

A surprise comes about when we build the matrix valued polynomials going 
along with the scalar ones alluded to above. 

We consider the differential operators with matrix coefficients and order not 
exceeding three. Notice that for generic values of a, (3 this gives only the iden- 
tity operator. The situation now is entirely different. We find four nontrivial 
operators, which along with identity, give a basis for this five dimensional space. 
The matrix valued polynomials P„ (x) are common eigenfunctions of two linearly 
independent third order operators. But there are other surprises in store. They 
also satisfy differential equations of orders one and two, namely, recalling that 
/? = a + 1, 



x-1 




+Pnix) 



1 



R_ 

2/3 



/3- 



2/3 



1 + 



2/3 



JL 

2/3 



1 




as well as 



, /O x(x-l)\ r-,/ f \ x + a \ „ , , . 



2 






^0 + n ( 2^+3 1 + 



1 
1 



Pn{x) 



where 



4(a + 1) V 
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